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Abstract 

Two new sets of models for describing compactly the Fourier Transform 
of Kohlrausch-Williams-Watts, both based on the adiabatic variation of 
parameters of a double Havriliak-Negami approximation along the whole 
interval of frequencies, are presented. One of them is relying, obviously, 
on the use of a well-behaved-pair of patches of the mentioned type of 
approximants, Ap 2 HN{uj). The other is obtained by altering the simple 
functions HN{uj) and making dissimilar the couple. They are proposed 
the guidelines of a new and systematic approach with extended Havriliak- 
Negami functions which is global, (non local), and of constant parameters. 
The latter at the cost of a more complicated dependency with the low 
frequencies than 1 + {iuJTHN)°‘■ 


Introduction 

To the extent that the object of study of soft matter and fluids has been passing 
from simple polar liquids to polymer, glasses and quasi-amorphous materials, the 
phenomenology of rheological, or dielectric, relaxations of physical systems has 
become increasingly complicated. And it is not only that several types of those 
relaxations are superposed along frequency space making difficult to distinguish 
among them, but that employed functional form evolves from an easy one as 
Debye [1], to other more complex as Havriliak-Negami |2j[3], > 

after experimenting with intermediate stages as Cole-Cole |1], > and Cole- 

Davidson |S], 

* tlazcala@yahoo.es 
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Simultaneously something similar happens in time description while we con¬ 
sider the different temporal scales implied, so certain habit to model, -imposed 
by the mentioned physical phenomena and the discriminatory ability of ex¬ 
perimental equipment-, shows a methodological exhaustion. In this sense any 
testing for the use of new relaxation functions, giving account of the novel ex¬ 
perimental records, is fully justified 017]. 

However, as has been quoted previously (liini nnuni nn usmi usi), the 
swapping from one functional space to the other still continues to be hard, 
often drawing upon, the researcher, efficient numerical methods to perform such 
devious change [T5J [T7J [TBIIT^ . 

Nevertheless it turns to be sometimes unsatisfactory this capability for blind 
calculations as it does not provide many times of a general view allowing for the 
interpretation and identification of these models whose available information 
is fragmentary or incomplete. In this sense a catalogue for formulae linking 
frequency and time realms mi, is a precious help while are appraised significant 
system parameters. Besides it is also an essential partner of the numerical 
analysis to bound errors and expose constructs [SIIISIIII], both characteristics 
of computer techniques. 

We already have shown as a set of Weibull distributions [22], exp 

0 < < 2), in Fourier space, admit a good approximate description by 

sums of Havriliak-Negami functions [231 123] [53] . Additionally in a precedent 
paper to this work | 2 S], it is established the local character of the approxima¬ 
tion and how, with slight variation of the parameters {ai, 2 , 71 , 2 , ti, 2 ) A} with 
frequency ui, Ap 2 HN can describe a perfect fit with the objective function, tpp. 
Such adiabatic behavior is commonly misunderstood as an argument against 
the approximation by means of basic relaxation functions as Havriliak-Negami. 
This fact it is best interpreted as the need for a wider family of relaxations with 
a known local portrayal. 

Thus in this job we will focus on taking advantage of such “local” information 
to build “global” functions that ameliorate the preceding approximation in the 
whole range of frequencies, [0, 00 ). Also the relative error of all proposed models 
in the present and previous report | 2 S| is depicted and tested against the real 
data obtained from Fourier integrals. 

1 Uncommon approaches 

1.1 The global two-term approximant: version with an 


atlas 

In short we have constructed two approximants of type |25| . 



( 1 ) 


(Ai -I- A 2 = 1), for two different overlapping intervals v G [0,500.0005] and 
V G [1.00,10^^], (or ly G [1.00,10^] if /3 > 1). Besides if we consider how the 


2 



relative error between moduli of approximant and function behaves as frequency 
varies, (be. it stabilizes at an almost constant value never greater than 0 . 2 % 
for high frequencies and 0 < /3 < 2), the upper bound of second interval can be 
extended without a big amount of error to an unlimited frequency. It means that 
there are two charts {Ap 2 ,iHN{oj),Qi) and (Ap 2 ,hHN{u;),nD with = 27rx 
(0,500.0005) and fl)) = 27rx (l,c») that reconstruct in an acceptable way the 
function in the whole interval ( 0 ,cx)), plus the value at '!/’/ 3 ( 0 ) = 1 as an 

imposed condition. 

Now, all that is needed to obtain a global solution is a way to merge both 
charts without overlaps, following a standard procedure. We will resort to an 
smooth, and monotonously increasing, function defined as 0 

{ 0 : uj < oji 

( 0 , 1 ) : oji < UJ < LOs , 

1 : ojs < UJ 

with uJi and Wg chosen arbitrarily. So starting from locally adjusted functions 
which are properly selected it is possible to write a suitable approximation to 
function ijjp, in the whole interval [ 0 ,oo), as: 

Ap 2 HN{uj) = Ap 2 .iHN{uj) x (1 - W,.g) + Ap 2 ,hHNiuj) x ( 2 ) 

With our particular choice uJi = 2 tt and uig = 47 r, (or uig = 2.27r if /3 > 1), 
we have finally laid down a global surrogate of Havriliak-Negami type for the 
Fourier Transform of Weibull function iljp{uj), whenever 0 < /I < 2. 

1.2 The global two-term approximant: version with adia¬ 
batic parameters 

1.2.1 The stretched case /3 < 1 

While the chart for low frequencies is a quite good approximation to function 
tppiuj) it is not however a perfect matching in the mentioned range. Thus, 
there still is room for improving the fit as the difference reaches a peak of 1.3% 
for some of the lower frequencies, uj < 10<5a;, at medium values of beta, i.e. 
P « 0.30, and working with an r 2 sampling, (Sui = 1/999.999). So the logical 
next step would be to add another Havriliak-Negami term to the approximant 
to fill the gap, although using some restrictions, (v. gr. over a ■ 7 products), 
to avoid proliferation of parameters and to hold the resemblance of a truncated 
series. The procedure works, however one should cope with other problems 
of the multi-parameter optimization as the non isolated loci of minima, the 
multiplicity of them, and the competition among coefficients for some regions, 
(which prevents a broad and proper allocation of them in the search space). It 
is then a good idea to analyze why a small gap befalls precisely at very low 
frequencies in order to correct it or design a strategy for the restrictions of new 
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Formulae 

/3A expE^— 

a,{i - ar] 

)^] dsC - /3)^ 

Parameters 





Constants 



Parameter 

1 A 

A 

0.886396 

1.55239 

Constants 





B 

0.186694 

ai 

0.606292 

0.282576 






di 

1.65185 

(^2 

1.46112 

6.08038 







-4.51252 

as 

-8.94199 

-32.0705 







3.81826 

a4 

13.2755 

47.4824 







-0.916873 

as 

-6.3231 

-22.507 






Corr. 

0.996672 

Corr. 

0.999970 

0.999913 




Table 1: Optimization parameters of the Eq. [g] approximant, {case j3 < 1): Formulas, (with 
e = 10“^®''), and their constants for ot 2 , «3 and 1 — A. 


incoming terms of a series, or an approximant. 

The first derivative of with respect the angular 

frequency is written as: 


lim = lim -i 


(3tPe 


*dt = —i 


* dt = —zr(^ + l), 


,1 


w—J-0+ dcJ a;—J-0+ jQ 

while the first derivative for the Havriliak-Negami function is expressed as: 


d 1 _ -7a(iT)“u;“-i 

dw (1 + (*Wt)“)'>' (1 + (itLiT)“)T'+^ 


0 

—ijT 

—i°‘oo 


a > 1 
a = 1 
ct < 1 


which implies, as —= — cos ^ — i sin ™ takes values in the third quadrant, 
that no lineal combination of two vectors 


dHN.iu:) , 

( -A- 

dui 


)u>=o and ( 


dHN2iu}) 
dw 




can give a vector parallel to —i when both ai ^2 < 1- Nor it is possible with 
only one a < 1. In such a situation, i.e. ai < 1 or 0:2 < 1, the magnitude of 
at least one modulus will be infinity due to factor what gives a reason, 

-the momentary diminishing of any \HNi^ 2 i^)\ in the vicinity of w « 0“*" is 
quicker than that of |'(/'(w)|-, for the underestimation of tjtp by Ap 2 HNa,-y^T,x in 
the range of low frequencies. 

All this confront us with the fact of finding alternatives to 0 : 1^2 < !> and 
under the conditions of module and direction of the derivative oi tpp these are 
simply ai _2 = 1, or ai = 1 and 02 > Ij (or the converse pair of indices). 
That contradicts the empirical finding for parameters ai ^2 we made with the r 2 
sampling in the range of medium to low frequencies which supports heavily the 
condition 01^2 < 1 when P < 0.60 and the conditions ai < 1 and 0:2 ~ 1 when 

practical production of function can be made following instructions found in 

ref. |26], lemma 1.10, p. 10. 
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<1 - \ + Ss = l co=(C,3 + 0a)} 

^r=Q 


Parameter 

Constants 


loglO "^1 


^>0 

0.0979577 

bi 

0.757724 

^2 

-0.228709 

Cl 

-0.142877 

Cl 

22.2049 

A 

-7.14011 

C 2 

0.649622 

C 2 

-8.86238 

4>2 

5.19663 

Corr. 

0.999836 


exp [(^)=E?=o e»(l-3) = ] 


Parameter 

Constants 

r2 

CO 

-0.796921 

Cl 

4.14764 

C2 

-2.98386 

C3 

-14.9068 

C4 

33.8185 

C5 

-26.9545 

C6 

7.67633 

Corr. 

0.995682 


exp[- M0^] 


Parameter 

Constants 

loglO t-3 

M 

2.74977 

d 

4.40586 

bo 

-9.15927x10“® 


0.559288 

b2 

0.159154 

h 

-0.641584 

Corr. 

0.999992 


Table 2: Optimization parameters of the Eq. approximant, (case /3 < 1 ): Formulas and 
their constants for ri, T 2 and 


j3 > 0.60. Consequently the double approximant will always hold this gap in 
the region of small frequencies unless of course we add some local modification 
to parameters {a, 7 , t, A}i _2 near w « O'*'. 

In summary there are three zones in the w-space where the coefficients 
{a, 7 , r, A}i _2 are similar in magnitude, -comparing equal symbols-, although 
with a different behaviour as functions of /?. And they change from one to 
another conduct when they are shifted across the intervals of frequency. This 
means the set of parameters depends on w, i.e. {a, 7 ,r, A}i ^2 = {a, 7 , t, A}i_ 2 (w), 
but they are of a very slow variation through the whole interval of frequencies 
[0, 00 ). They are adiabatic coefficients of the approximant Ap 2 HNa^^^T,\i^)- 

Therefore we are looking for a function reproducing the main traits of 
namely: (^)<.p=o II * |#(w)| 0 ( 1 ), and -i/’/ 3 (w) « 0 {Ap 2 HN{uj)) 

locally Vw S [0,oo) with adiabatic coefficients as parameters in the fit. So with 
that in mind, and respecting conditions (a • 7 ) 1,2 = P for /3 < 1, our candidate 
will take the form: 


1p/3{uj) = AMi^2HN{uj) = 


A 

(1 + iWTi)^ 


1--A 

(1 + Ml(uj){iUJT2)°‘^)^ 


(3) 


with the mollifier of the Havriliak-Negami function, verifying the fol¬ 

lowing boundary conditions 


lim Mi{uj) = 1 and Mi{uj) « (wts)"^, j/ w « 0, 

UJ —>-00 

where a 2 + 0^3 > 1- 

Undoubtedly to determine the mollifier it is not an easy task and surely its 
expression as a series could be as difficult of obtaining as the one of ipp, and yet 
the conditions we have imposed on may allow an easy estimate of it. So 
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Figure 1: Optimization parameters of modified Havriliak-Negami approximant that is shown 
in Eq. They, 0 : 2 , ti, t 2 and 1 — A, follow strictly Eq. [^while T 3 and 0 : 3 , ( 0:3 in text), are as 
in Eq. the latter a simple estimation for the mollifier. 0 < 1. Solid lines are mathematical 
adjustments for parameters, which are given in tables^ and 


we will put forward as an estimator of mollifierA^;(a;) the function 

Mi{u}) = [- arctan((a;r3)“^)]^, (4) 

TT 

setting TV = 3 as an appropriate average after a timely optimization for some 
values of /3. 

1.2.2 The squeezed case P > 1 

As ai > 1 and 02 < 1 for 2 > f3 > 1, when the adjustment to r 2 sam¬ 
pling is done, it is arguable that a lineal combination of both direction vectors, 
( dgjVi^ 2 (<^) at w = 0 can be parallel to —i, since they take values in 

opposite quadrants {i.e. second and third ones). Nevertheless the magnitude of 
the derivatives, (zero and infinity respectively), avoids such result, and again the 
admissible options for 012 are those of the case /3 < 1 which as we pointed out 
contradict the numerical findings. We must proceed again with some modifica¬ 
tions of the approximant, however it is not possible now to extend the model 
of equation to this situation. Firstly the logarithmic second derivatives of 
logio|' 0 / 3 | are qualitatively different for /? < 1 and /? > 1 cases. (See figures 
1 and 2 in [IS]). Second only the 02 • 72 = /? condition holds while the other 
becomes oi • 71 ~ 2/3 for the asymptotic behaviour of tails. Besides in the low 
to medium frequencies range the products (a • 7 ) 1,2 present a nonlinear trend 
greater in many cases than the corresponding linear condition of high frequency. 
(See figure 9 in |25| 1. In consequence we introduce a peculiar model with two 
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Figure 2: For /3 > 1 an approximation to is done using the modified approximant 
set in Eq. At left panel parameters: tq, "ti, T 2 , and A. At right panel: ap, 

Q 2 , as and 71 . Except for rs and Q 3 , ( 0:3 in text), all of them are defined in Eq. The 

former are established in Eq. which describes an estimator of the current mollifier. As in 
previous graphics the parameters are adjusted with suitable mathematical expressions, (solid 
lines), given in tables [^and[^ 


characteristic times in the first Havriliak-Negami relaxation which jointly with 
the two terms structure of the whole approximant will retain the mentioned 
characteristics and restrictions for tails, with the sole exception of oq • 71 which 
is to be determined by fitting. The latter in light of the numerical results ap¬ 
pears a spurious property or at least a virtual one conceived to explain the 
sudden change of curvature in a small interval of frequencies, v 

Thus when /? > 1 the global formula for the approximant it reads: 


'0/3(a;) = AMg^2HN{(jj) = 


_A_ 

(1 -I- iojTi + 


with the boundary conditions for the mollifier as 


1-A 


(1 -I- 7Wg(a;)(ja;r2)“2)<^2 

(5) 


lim A4g(uj) = 1 and Mg{uj) ~ 0 ( 1 ), if uj 0, 
oj —^00 

where oq > 1 and 02 ^ 1- 

Again an exact expression for the mollifier is out of scope of present study 
and we settle for an estimator like 

Mg(ijj) = [^^{1 + (— 1)^ arctan(a;T3)“^}]”, (6) 

V2 

and with an ad hoc choice of n = 3. With such estimator we will obtain a good 
approximation to \'4ip\ which deviates slightly in a neighborhood of w ~ 27r, (the 
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zone where the maximum of curvature happens in logarithmic scale), although 
it describes fairly the body of function and quite well the trend and values of 
tail. (See figure 2 in [25]h 

1.2.3 Analysing trends of parameter graphics 

Finally in figures and we display the parameters of estimated functions 
AM.i^ 2 HN{uj) and AM.g^ 2 HN{uj) and their adjustments as functions of /3 (re¬ 
spectively /3 < 1 and /3 > 1). In this occasion the interval of frequencies is up 
to V = 10 ^^, and up to z/ = 10 ® depending upon choice of /3, and the sampling 
of frequencies is what we called logarithmically homogeneous, i.e. Tsi. In both 
cases all the curves have a break, or turnaround, more or less evident according 
to each one. This occurs for each curve, -within same j3 case-, in the same point 
(/3 « 0.80 and /? « 1.80). Also during the course of several optimizations such 
points have changed marginally and the breaks have increased or diminished 
their sharpness according as we changed the average exponent (A, n) of equa¬ 
tions]^ and data weights or samplings (r 2 , Vsi). So in conclusion we interpret 
that such abnormalities are a result of the shape of estimators. 

Stretched instance, /3 < 1 

The best option in case /3 < 1 is to weigh, -while using xmgrace to get a fit 
m-, the tails with option l/T^ to soften the jump and obtain an even ad¬ 
justment all the way in the interval of frequencies. The value of N also could 
be lowered but the price to pay is an increasing error for all the matching be¬ 
tween both functions, (approximant and ipp), around values of /3 G (0.1,0.3) and 
/3 € (0.7, 0.8). On the other hand the ability of the new function AM.i^ 2 HN{uj) 
for describing the effect that slow variation parameters {q;, 7 , r}i _2 would have 
in the original Havriliak-Negami functions fully justifies the introduction of mol- 
lifier Unfortunately the expression of its estimator does not seems good 

enough in the vicinity of w « 0 since the results do not fulfill the required condi¬ 
tion 02 -l- N&s > 1 at all when /3 —> 0“*". This is a consequence of having frozen 
the exponent at A = 3, we should increase its value till infinity to compensate 
the empirical trends of 02 and 03 to be zero when /3 —)■ 0^. However the first 
term of AAii^ 2 HN{uj), an almost residual one since A « 0 for /3 < 0.8, seems to 
balance numerically this mathematical unsuitability of the second term of the 
approximant in the description of ipp{uj « 0). And that is possible since there 
is no conflict in accounting for a slow diminishing I'lppl in the neighborhood of 
w « 0, (/3 < 0.10), using a fast decaying Havriliak-Negami type function, {i.e. 
of large r), with the sampling step that we used. For such small values of (3 
a neighborhood of zero where I'tppl ^ 1 is so elusive that a frequency step of 
Siy = 10 “® is too large for considering a description of the modulus gradual 
decay. (In tables and we wrote the mathematical expressions for the six 
parameters of A(tLi) as curves depending of variable /3). 


Formulae 

exp [{Zs^o 0 . 3(0 

- 1)^ } exp(-M^)] 

B + exp[ —— I)"- 

lEs = l f>s(/3 - 1) = 

C0^ expE^^ 

- 1)”] 

Parameters 

QQ ~ 2 

0.880214 

02 — 1 

4.1814 

Parameter 

0:3 

1.45629 

Parameter 

Constants 

c 

Cl 

71 

0.730891 

3.50858 

Constants 

M 

Constants 

M 

ao 

1.70465 

0.0522961 

B 

2.0503 

C2 

-34.1656 

ai 

0.976489 

-0.770667 

bi 

-9.16895 

C3 

153.848 

a2 

2.74286 

2.78766 

b2 

30.8241 

C4 

-373.845 

0.3 

-2.8097 

-2.45835 

h 

-41.314 

C5 

507.468 

a4 

1.36045 

-0.148637 

b4. 

22.1584 

C6 

-359.883 

Corr. 

0.963441 

0.98941 7 

Corr. 

0.998267 

C7 

102.995 


Corr. 0.999860 


Table 3: Optimization parameters of the Eq. approximant, {case /3 > 1): Formulas and 
their constants for aq, 021 “3 and 71 . 



BJtp [(-g 


- 1)"] 


exp [ - M(f3 - 1)«1 

93(i> - i)“ 

Parameters 

Constants 

P 

do 

TO 

= 11.5 

1.40526 

Tl 

= 11.5 

1.91075 

T2 

= 4.4 

-0.0241828 

T3 

= 9.8 

1.36359 

Parameter 

Constants 

M 

A 

5.82578 

di 

-38.3857 

-27.9873 

-4.7472 

-30.3907 

91 

273.583 

d2 

503.134 

484.97 

10.0218 

380.42 

92 

-407.192 

dz 

-3368.51 

-2810.17 

-31.8961 

-1955.88 

93 

-3.32942 

dzi 

8385.44 

7209.62 

44.4373 

3128.49 

94 

195.919 

ds 

-9737.49 

-7533.20 

-28.0215 

-2939.76 

Corr. 

0.998083 

Corr. 

0.999708 

0.999695 

0.999214 

0.999346 




Table 4: Optimization parameters of the Eq. approximant, (case /3 > 1): Formulas and 
their constants for tq, ti, T 2 , T 3 and A. See figure]^ 


Squeezed instance, /3 > 1 

Case /3 > 1 is instead more difficult to adjust in the whole interval of frequencies 
since no additional weight is possible to use. The kink of log^Q near v ^ I 
claims for a body not overlooked which would be the case if tails were given more 
importance by weighing them as in previous procedure. Besides, the mollifier 
of Havriliak-Negami function is not enough elaborated and as a consequence 
appears a bifurcation for each curve of parameters, corresponding the lower 
branch to the best adjustment to data. Nevertheless if the latter is employed for 
describing the curves, an abrupt change in trend for them is evident and makes 
more difficult the handled mathematical expressions in parameter adjustment. 
We show here only the upper branch of all curves, this leads to a smooth and 
nice interpolation line for each parameter as seen in tables and 

Although we started with a nine parameters ansatz for the approximant 
AM.g^ 2 HN{u}) is clear from the graphs, (right panel of figure]^, that uq and 
a 2 are almost constants. Now the written requirements over them are amply 
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fulfilled. Only there is a small disagree of order 10“"^ from condition a 2 {P) — 1 
for some values of (3. This is entirely due to competition among parameters 
and subsequent numerical errors. Meanwhile ao{/3) « 2 for all betas, and any 
of them differ from this number less than 1.5% for (3 > I and only with some 
significance for beta 1.00 and 1.02, and for (3 > 1.80. Thus with a slight setting of 
A4g(uj) has to be possible to write AA4g^2HN{uj) as a seven parameter function 
which is more economic computationally. (All the adjustments to this set of 
parameters are given in tables and . 

1.3 The role of the ’aide-de-camp’ in the modified approx- 
imant 

Apart from already explained conditions in the onset of frequencies which makes 
a Cole-Davidson relaxation suitable to describe the boundary condition of 
it is obvious that in the interval f3 G (0.80,1.00] the first term of approximant 
described in Eq. plays an important role in the approximation since A is not 
at all negligible. However for the interval f3 G [0.02,0.80] the story is quite 
different, almost all its contribution is forced by theoretical considerations as 
now the share coefficient is really small. To extend this situation and make an 
adjustment with a one-term approximant in the whole interval f3 G (0,1), we 
should prepare a more flexible second term of Havriliak-Negami type in Eq. 
And with this goal in mind we establish the exponent N of as a new 

parameter of the optimization. Namely we shall do the following setting: 

V'/3(w) = AMi^iHN{u}) = -- -- -—— 

(1 -I- Ml{uj){iUJT2)°‘^)‘‘2 

with an estimator to similar to that of Eq. [^though now N = N{(3) is 

non constant. 

The results, {i.e. the parametric curves of /3), are shown in figure]^ there 
we note two important features about the behaviour of parameters a 2 and da 
and the shape peculiarities of N{j3). The first characteristic, in the interval (3 G 
(0.80,1.00], is that we do not recover the functional form of a Debye relaxation 
as ,5 —>■ 1. For such a requirement it should happen at least JV 0 and a 2 —t 1 
as an strong condition, or 02 + Nas ~ 1 as a weaker one, in that limit. Neither 
the strong nor weak conditions are fulfilled by the parameters as can be seen in 
right and left panels of figure [^ 

In light of the share coefficient behaviour (A) remains an important question: 
if the auxiliary term dominated by it in the modified approximant is really 
necessary. (See left panel in figure Q. Or if instead it is only needed to ’unfreeze’ 
the exponent N in the estimator A4i{uj) of the mollifier, (see Eq. |^, to adjust 
properly with only one term: the mollified Havriliak-Negami function. 

The second flaw is patent when we realize that it is not possible, in the 
interval (3 G (0.02,0.10), to hold the condition 0:2 + > 1 when 02 —t 0+ 

and 0:3 —>■ O'*' since N is finite and decreasing as /3 —>■ 0. These trends of 
alpha parameters are attested, jointly with the N one, and depicted in figure 
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Figure 3: Fitting parameters of Eq. [^as functions of variable 0 <1. Left panel: character¬ 
istic times T 2 , ra and exponent N{0), (see Eq. |^. Right panel: frequency exponents Q 2 and 
“ 3 , (03 in Eq. The dotted lines are just guides for the eye. 


again. In conclusion the modified relaxation of Havriliak-Negami fails in 
the adjustment at both ends of jS interval and it is not hard to imagine the 
difficulties it has to describe an environment of \ip 0 \ ~ 1, (*-e. w « 0), with 
a poor sampling of very low frequencies, (as is the case of ours for so small 
values of beta). The tails obviously, in such a situation, lead the adjustment 
and the mentioned requirements about the behaviour of near w « 0 should 
be imposed externally. At this event the best option to save both flaws is to 
maintain the optimization with a two-terms approximant like that of Eq. 

At last a further reason to keep the expression of Eq. is to have a formal 
similarity with Eq. for the case /? > 1 which makes more manageable the 
treatment and analysis of the problem in the whole interval /3 G (0, 2). 

2 Comparison and discussion: Suitability of for¬ 
mulae 

In light of these circumstances we will consider the frequency-averaged relative 
error, (among data and test functions), as an indicator of reconstruction 

capability for any of the proposed Havriliak-Negami approximants. As it has 
been patent till now most of the present discussion here refers to the suitability 
of pairs proposed to describe the modulus of data We must be aware that 
aside from the results here discussed some additional tuning of phase should 
be sought. Different one with each model for approximation we use. Even so, 
without all the benefits of the phase, an accurate adjustment between data and 


11 









Figure 4: Frequency-averaged relative error, , for data ipjs = 

obtained with Matheniatica'^^ and function values of models AAC (circle), Wutt (dot- 
dashed), AtAdj (plusses), Atint (squares), HTIt (light solid), GlbAdj (dark solid) and Glblnt 
(diamonds). Eight frequency windows, whose details are written in the text, and two types 
of samplings -linear from (a) to (d) and logarithmic from (e) to (h)- are shown for interval 
0 < /3 < 1. 


approximants makes this methodology of multiple Havriliak-Negami summands 
useful to determine form parameters in dielectric spectroscopy experiments, or 
in systematic search of them by means of genetics algorithms. 

Previously, a frequency-averaged relative error for the moduli of functions 
was defined as: 


Such error function is depicted in graphics and There it was calcu¬ 
lated for the different models described in a previous paper, |25], and in the 
present one, and also for other two models found in literature (see references 
HHHD]). In particular, the errors for a Havriliak-Negami function, the three 
models described above (double H-N, atlas of aproximants and modified H-N), 
two variants of the last ones with parameters calculated via formulae given in 
tables 1 to 4 in [25], and the numerical solution obtained from the C code in 
reference m are given, when /3 < 1, in figure Same models with /3 > 1, 
except IHN of AAC (Ref. [20])) are portrayed in figure]^ 

The frequency windows studied are: (a) [uJm,^x )= [0,1)) (b) [1,10), 
(c) [10,100), (d) [100, 500] in the upper row of both figures. They followed a 
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linear sampling with 5v = 999 ^ 999 , conversely the lower rows were logarithmic, 
homogeneous in each decade in the way we already explained. Their intervals of 
frequencies are: (e) [ 100 , 10 ^), (f) [ 10 ^, 10 ®), for both graphs, but (g) [ 10 ®, 10 ®), 
(h) [10®, 10^®] with /3 < 1, and (g’) [10®, 10^] when /3 > 1. 

For a quick sight inside the plots in and we have tagged the models 
already explained. We remind that obtained from the 

direct calculation of Fourier integral and is the same reference for all errors 
{Srei) calculated with different test models which are now listed as: AAC, 
the Havriliak-Negami approximation cited in ref. [5D]. Wutt, the C library of 
reference m which employs the power series for low and high frequencies and 
an effective numerical method for the intermediary frequencies in the interval 
/3 € [0.1, 2.0]. AtAdj is the label assigned to the model of equation]^ and the 
same is true for the symbol Atint. The distinction is that while the parameters 
{ 0 : 1 , 2 , 71 , 2 , 7 ‘i, 2 , A} in the first case are calculated following the formulas in tables 
1 to 8 of | 2 S|i in the last one are obtained from the points directly obtained of 
error minimization and depicted in graphics 5, 6 , 7 and 8 of |25| . For the 
sake of clarity we have repeated the latter results showing separately each part 
which follows the eq. for low or high frequencies (head and tail functions 
of equation [^. It allows to appreciate where in the frequency interval the 
individual approximant diverges from the data, and which one is exactly its 
contribution to the atlas of approximants. The transition from a plain downhill 
to a potential tail it is then quite clear. This is called HTIt. Besides all 

the three previous models refer to the range [ 0 . 02 , 2 . 00 ] of shape parameter /?. 

The last assertion it is also true for labels GlbAdj and Glblnt though two 
different formulas and their respective implementations are employed, the equa¬ 
tions [^ and [^ for j3 <1 and the equations [S] and for (3 > 1. Again the first tag 
refers to the adjusted parameters (see tables [T]^ and and the second to 
the original points as depicted in figures [T] and^] 

2.1 Models response in the stretched case 

In figure [^ A AG, the approximation with only one HN function, shows the 
biggest of all errors for the models here presented and the interval 0 < /3 < 1 . 
The best result is of course that of Wutt which combine analytics and numerics. 
Meanwhile the atlas described by equation]^ {i.e. Atint), works quite well even 
in the interval of very low frequencies where we demonstrated that one of the 
functions of the approximant should be a Cole-Davidson relaxation or a modified 
version of it, -this depends on if /? is less than or greater than unit. (See Eqs. 
and|^. Besides, it holds good terms in the medium range thanks to the change 
of describer function, (i.e. from head, Ap 2 ,iHN, to tail, Ap 2 ,hHN-, see figure 
Impanel (b)). The performance of this swap is even better for high frequencies 
because the matching with data exceeds expectations and the approximation has 
not required any restriction in the product (q:- 7 )i_ 2 - This reinforces our previous 
conclusion |24| . and points to the true nature of as a sum of a Havriliak- 
Negami pair with almost constant coefficients which only change significatively 
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as 1 > w —>■ 0 pS]. Such important transition is highlighted in figures]^ andby 
the discrepancy between model Atint (squares) and HTIt (light solid line) and 
reminds us the need for at least two charts (one for z/ < 1 the other for z/ > 1) 
in the description of the whole function ipp. Usually this approach is not taken 
in consideration in the literature since only one set of parameters {a, 7 ,T, A} is 
employed to match the data, or if considered is misinterpreted due to the usage 
of a time scale factor in the stretched exponential {i.e. exp—{t/TKww)^) |29| . 
See panels (a) to (d) in figure 

It is worth to note how the Havriliak-Negami approach A AC starts to work 
better than the double approximant HTIt in the range of medium-to-high fre¬ 
quencies, according as 0.8 < /3 1. It sounds logical since /3 —>■ 1 implies A —>■ 0, 
and this final value is a pathology for the double sum of HN functions quite dif¬ 
ficult to treat numerically. Such problem is not evident using the atlas Atint 
since the approximant of tail Ap 2 ,hHN takes the control over Ap 2 ,iHN in such 
frequency interval. The former function besides, at very high z^’s and with /? 
near 1.00, shows similar errors to the results of numerical-analytic method Wutt. 
Now the major problem for both of them will be the numerical oscillations of 
reference data. See panels (g) and (h) in figure]^ 

The differences of error between model AtAdj and Atint present clearly two 
regions. One in the low frequencies zone {v < 1) the other in the high frequency 
one (z^ > 1), as it is usual coinciding with the onset of potential behaviour for 
tails. In the first case it is the lack of ability of the double approximant with 
constant parameters to approach data, what gets closer both models. This is 
shown in figure and is more conspicuous when ,5 ^ 0.40. In the second case 
where both models split apart more than one order of magnitude the reason for 
this is more subtle because the correspondence between and is tighter 
as both functions follow a quite similar potential decaying. (See graphpanels 
(b) to (h), and figure 9 in [2S]). What makes such difference between Atint and 
AtAdj is an extra error provided for each of the parameters ai, 2 i 7 i, 2 i ti ,2 and 
A. An additional contribution which is consequence of the fitting of parametric 
curves to optimized points. Then it would be desirable to reduce those inputs 
binding the parameters to relationships as the already mentioned (a • 7 ) 1,2 ~ P- 
Nevertheless, we feel that some further work should be done to link the conduct 
of and A to the coefficients of the analytical series for and so in the 

absence of them we have presented the models drawn from eq. free of any 
external constraints. 

As a novelty we introduced here two global models to simulate the data, 
namely Gbllnt (optimized points) and GblAdj (adjusted parameters). Watch¬ 
ing them carefully in the various intervals of frequency one can realize how the 
performances are similar to the models of Havriliak-Negami double approxi- 
mants, Atint and AtAdj, respectively. Also is possible to observe how for the 
low frequency range, n < 1, the model Gbllnt outperforms to Atint, although 
this one later improves and is usually better for the high frequencies in accuracy. 
(See in panels (a) to (h) of figure]^ the square and diamond curves). Also for 
the global model, the aggregated error of several independent parameters spoils 
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Figure 5: Relative error, for data obtained with 

Mathematical and test values ^ of models: Wutt (dot-dashed), AtAdj (plusses), Atint 
(squares), HTIt (light solid), GlbAdj (dark solid) and Glblnt (diamonds). Seven frequency 
windows and two types of samplings -linear from (a) to (d) and logarithmic from (e) to (g)- 
are shown. Maximum frequency in panel (g) is reduced to Ux = 10^ due to glitches and 
numerical noise of data beyond this point. 1 < /3 < 2. 
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the result making alike the error curves of cases AtAdj and GblAdj, (plusses 
and dark solid lines in graphics). Therefore the drastic displacement of rela¬ 
tive error lines towards a lesser precision for models with adjusted parameters 
7 i,2j Ti,2y '^}(/ 3 ) points to the need for strict relationships among them and 
good descriptions of functional dependence with the shape parameter / 3 . 

However it is important to emphasize how proposed estimators of mollifiers, 
in Eqs. and [^, are subject to many “ad hoc” restrictions, mainly 
deduced of data traits and information obtained from the behavioral changes 
of 01^2, 7 i,2, ti,2, a curves while changing the regime of frequencies from low to 
high. The most significative restraint here is the fact the AAi^g functions are 
set as real ones when they should value in the complex field. As we have seen 
before, in Eqs. and near w « 0 ^ the dominant demeanor of Mi^g{(jj) is 
determined without consider further modifications to the phase since the 
frequency factor 0;“^“'"* extinguishes such contribution quickly and only remains 
the one of extended Cole-Davidson term. Nevertheless the role the phase plays 
is over the entire interval of frequencies and, although to our present purpose 
of describing the modulus accurately it is not crucial this bias in the argument 
of the approximation AMi g 2HN, it is important to point out the need of a 
mollifier in complex series for proper description oi ipp- 

It seems promising, from an analytical point of view, that a strategy to sum 
up a quite difficult series in the neighborhood of w « O’*" comes from the help of 
an extended Havriliak-Negami pair. Perhaps could be interesting to decide the 
mollifier’s functional form with the aid of series, integrals or equations which 
determine ipp. Mainly when w —> 0 , for /3 < 1 , or when w —> 00 for /3 > 1 , the 
most difficult cases for power series involved Emails]. 

Obviously, in the light of the problems we face when using unsuitable func¬ 
tions as estimators of J\ 4 i^g{uj), it is suggested that a similar ’loss’ of phase could 
happen in the atlas approximation of eq. to data. And consequently foresee 
a mild shift in the argument of whole function AP2HN{oj) at high frequencies. 
This displacement is due to the way the tail functions Ap2,hHN{ijj) are deter¬ 
mined. Data are pruned in a logarithmic pace and the important information 
at low frequencies, -the plateau-, is removed when optimizing tail parameters, 
something is not made in the case of head functions, Ap2,iHN{uj). All this does 
not impact very much on the approach to the modulus, as we saw in graphics of 
figures 1 ^ and but suggests a share parameter fully complex, i.e. A S C — K. 
Unfortunately it would add a new degree of freedom and would overshadow the 
discussion about modulus characteristics in absence of a thorough treatment of 
the data phase. 

Apart its mentioned inability to describe the very low frequency range, the 
double Havriliak-Negami approximant, Ap2HNa^j^T,\i^), will not present such 
difficulties while describing the argument, -much less the modulus-, of data ipp. 

2.2 Models response in the squeezed case 

The last stay in this discussion is the figure that shows relative errors of the 
six previous models for shape parameter interval 1 < /? < 2 . As predicted by 
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first and second logarithmic derivative, (see frames a) and b) of fig. 2 in |25|'). 
there is a sudden change of behaviour in the slope of |'!/'/3| from flatness to a 
potential decline in a relatively small interval of frequencies. (See frame c) 
of same graphics in |25p. This is a much more sharper and distinct transition 
than in case /3 < 1, which forces the existence of a different set of constraints for 
products a • 7 in the double Hav.-Neg. approximation, as clearly shows figure 9 
in m- All this oblige to abandon quickly the HTIt model, (light solid line in 
panels (a) to (d)), in favor of Atint because the latter holds itself quite close 
to the potential tail and shares same description of with the former at very 
low frequencies. (See squares inside panels (a), and (b) to (g) in figurej^. 

Again, as with /3 < 1, a big distance in terms of relative error separates Atint 
and AtAdj, and as before this gap is attributed to a collective error subscribed by 
each individual parameter, whenever every uncertainty is caused by obtaining 
the appropriate parameter from a pertinent fitting function along all values 
of /3. However for the models Glblnt and GlbAdj such distance doesn’t exist 
at medium frequencies and thereafter, i.e. u > 1. (See diamonds and dark 
solid line in graphics of figure . It seems that model Glblnt it is not able to 
keep track of data tail so close as Atint does. Surely the ’bi-chronicity’ or the 
mollifier in functional form of Eq. [^should be revisited to give a proper account 
of directional twist of data near v ^ 1, and thus to diminish the error below the 
collective contribution of parameters. As in fact it happens at low frequencies, 
(see figure]^). Nevertheless, as far we know, this is one of the few attempts to 
describe globally the Fourier transform for /3 > 1 with an analytical model 
albeit approximate, so each piece of formula AAig^ 2 lIN{uj) has great value for 
future mathematical inquiries. 

Finally we should point out how, at very large frequencies, the potential 
description overcomes some numerical difficulties experienced by Wutt model 
for large /3 values, i.e. /? > 1.60. (See figures]^ and[^). 


3 Conclusions 

The present work is devoted to a compact description of the Fourier Transform 
of the Kohlrausch relaxation. As any reconstruction of this function as from 
spectral data should heavily depend on the information of frequencies near zero 
since tails will be surely corrupted by noise, an extra effort has to be made to 
manufacture a global function mimicking all aspects of this transform from low 
to high frequencies. Thus, two new sets of such models are proposed, and a 
detailed discussion on their errors compared to numerical control calculations 
generated directly by evaluating the Fourier integrals is made. 

We found how the approximation with a double H-N function always under¬ 
estimates ipi^ in the low frequency range, {v < 0.1). And although it is a small 
difference forces us to change those values of parameters {(a, 7 , t)i^ 2; A} already 
obtained in the range of medium frequencies. 

This repeats again in the transition from medium to high, or very high, 
frequencies. Nevertheless the double Havriliak-Negami approximation is close 
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enough to the original function as to describe it along a wide range of frequencies 
before the variation can be noticed. Moreover the parameters should not be 
regarded as varying, if the interval where the approximation is performed only 
comprises very high frequencies. 

Thus, due to the slow variation of parameters of the approximant with the 
frequency, {adiabatic parameters), instead of a global function with w—dependent 
{(a, 7 , r)i_ 2 , A}, we employed different charts of double Havriliak-Negami sums 
to describe locally -ippiuj). We found that using two charts is a good approxi¬ 
mation, enough to establish an atlas, however the differences at low frequencies 
still persist with such number of maps. The inclusion of a third one in the 
neighborhood of zero, i.e. v ~ 0, should be convenient, though the existence 
of an exact analytical series of powers in terms of Cole-Cole relaxations ^1E5]. 
{i.e. 1/(1 -I- {iu})^)), for ipp suggests the radius of such chart will depend on (3 
|29| . What it makes difficult and cumbersome working with an atlas of three 
maps. 

The question whether it is possible to sum up the series of Cole-Cole terms 
at iz = 0, or it is possible to write the atlas of Havriliak-Negami charts just in 
a global way, seems to have a positive answer. We presented two ansatze for 
extending the double Havriliak-Negami approximation, - that has proved to be 
successful locally -, which describe along several decades in frequency, and with 
enough functional proximity, the data of \ipp\. By means of rough estimates 
of the mollifiers of these new testing relaxations, we have found a good agree¬ 
ment with data moduli signaling a path for a future close approximation in the 
complex domain to the series, or numerical integrals, of the Fourier Transform 

V'/s- 
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